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Abstract 

Background: In recent years, there have been a growing number of studies on spatial inequalities in health 
covering a variety of scales, from small areas to metropolitan areas or regions, and for various health outcomes. 
However, few investigations have compared health status between cities with a view to gaining a better 
understanding of the relationships between such inequalities and the social, economic and physical characteristics. 
This paper focuses on disparities in respiratory health among the 55 largest French cities. The aim is to explore the 
relationships between inter-urban health patterns, city characteristics and regional context, and to determine how 
far a city's health status relates to the features observed on different geographical scales. 

Methods: We used health data describing hospitalizations for Chronic Obstructive Pulmonary Disease (COPD) as a 
proxy for respiratory health, and the total number of hospitalizations (overall) as a proxy for general health. This last 
indicator was used as a benchmark. A large set of indicators relating to socioeconomic, physical and amenity 
aspects of the cities (urban units) was also constructed. Data were analyzed using linear correlations and multiple 
linear regression models. 

Results: The results suggest that socioeconomic characteristics are major discriminators for inequalities in 
respiratory health status among urban units. Indeed, once combined to socioeconomic characteristics, only a 
climate indicator remained significant among the physical indicators. It appeared that the pollution indicators 
which were significantly correlated with COPD hospitalization rates loosed significance when associated to the 
socio-economic indicators in a multiple regression. The analysis showed that among the socio-economic indicators, 
an employment indicator derived at the regional scale, and two indicators reflecting the unequal intra-urban spatial 
distribution of population according to their education, were the most efficient to describe differences in the 
respiratory health status of urban units. 

Conclusion: In order to design effective urban policies, it is essential to gain a better understanding of the 
differences among cities in their entirety, rather than solely differences across small urban areas or individuals. 

Keywords: Respiratory health, Chronic obstructive pulmonary disease, Cities, Socioeconomic characteristics, 
Air pollution, Physical characteristics, Regional context, Intra-urban organization, Scale, France 
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Resume 

Contexte: Un nombre croissant d'etudes sur les inegalites spatiales de sante, conduites a des echelons tres 
differents, des quartiers aux aires metropolitaines ou aux regions a ete publie au cours des annees recentes. Ces 
etudes couvrent des pathologies variees. Peu d'etudes ont cependant adopte une approche comparative et ont eu 
pour ambition d'expliquer les differences interurbaines de sante en mobilisant a la fois des indicateurs sociaux, 
economiques et physiques. Cet article propose une etude de cas sur les disparites interurbaines de sante 
respiratoires parmi les 55 plus grandes villes franchises. L'objectif ici est d'explorer les relations entre les 
differenciations interurbaines de sante et les caracteristiques des villes, en tenant compte de ces caracteristiques a 
differents echelons (ville, region, quartier). 

Methodes: Nous avons mobilise une information sur les hospitalisations pour Broncho Pneumopathie Chronique 
Obstructive (BPCO) pour approcher la situation respiratoire des villes et sur I'ensemble des hospitalisations pour 
approcher la situation de sante generale de ces memes villes. Un grand nombre d'indicateurs socio-economiques, 
d'indicateurs decrivant I'environnement physique et les equipements de soins ont ete construits pour decrire les 
villes. L'information a ete analysee en utilisant des regressions lineaires simples et multiples. 

Resultats: Les resultats suggerent que les caracteristiques socio-economiques sont des facteurs majeurs 
de differenciation entre les villes du point de vue des situations de sante. Combinees aux indicateurs 
socio-economiques, seul I'indicateur des temperatures de janvier reste significatif. Les indicateurs de pollution, 
qui sont significativement associes aux taux d'hospitalisation pour BPCO, perdent de leur pouvoir explicatif lorsqu'ils 
sont associes aux indicateurs socio-economiques dans la regression multiple. Les resultats montrent que parmi les 
indicateurs economiques, le niveau de chomage regional s'associe a des indicateurs refletant I'inegale repartition 
des populations, selon leur niveau d'education, a I'interieur de la ville pour decrire de maniere efficace les 
differences interurbaines de sante respiratoire. 

Conclusion: Afin de mieux reflechir aux politiques de sante, il parait important de s'interesser aussi aux aspects 
collectifs des differences interurbaines de sante et pas uniquement aux differences entre individus. 



Introduction 

Epidemiological studies seek to identify and isolate factors 
that explain health inequalities in time and space, or 
among individuals. This study focuses on the spatial di- 
mension of health inequalities. There are different epis- 
temological ways of studying the relationship between 
health and space. One approach is to observe the individ- 
ual level in order to assess how far the social, economic 
and ecological characteristics of a place influence the 
health of its residents. Researchers adopting this kind of 
approach often use multilevel models. A second approach 
focuses on the health inequalities among different places 
and the resulting spatial patterns. The aim is then to iden- 
tify the relationships between health patterns and the 
spatial structure of different social, economic and environ- 
mental features. 

The research work presented in this paper is based on 
the second approach. It is thus an ecological study, and 
examines differences in respiratory health across cities. 
The aim is to explore the relationships between the re- 
spiratory health situation of different cities and their 
medical amenities, socio-economical and physical cha- 
racteristics. Intra-urban differences are also considered 
among the characteristics that can contribute to the over- 
all respiratory health status of a city. Various geographical 



scales are used in order to derive information on these 
characteristics: the regional scale, defining the context in 
which a city is located and that will be referred to as the 
regional context; the city scale; and finally the intra-urban 
scale, used to characterize the intra-urban organization of 
a city. 

Our investigation does not focus on the individuals a 
themselves. In this study cities are the objects of our study. 
They are not simply considered as aggregates of individuals 
but regarded as functionally coherent spatial entities, in the 
spirit of "cities as systems within system of cities" [1]. In- 
deed, as pointed out by Cummins et al. [2] we need to 
understand the complex and interdependent factors associ- 
ated with spatial health variations in order to design effect- 
ive "contextually sensitive" policy interventions to improve 
public health in cities. Our choice of focusing on the city 
level is one way of achieving this goal. 

The paper focuses on Chronic Obstructive Pulmonary 
Disease (COPD) which is a major cause of chronic mor- 
bidity (44 millions of patients diagnosed worldwide in 
2006) and mortality worldwide [3]. According to WHO 
estimates for 2030 [3], COPD is predicted to become the 
third leading cause of death in the world. Among respira- 
tory diseases, COPD is the most highly correlated with 
air pollution and its relation to urban pollutants is well 



Aschan-Leygonie et al. International Journal of Health Geographies 2013, 12:19 
http://www.ij-healthgeographics.eom/content/1 2/1/1 9 



Page 3 of 15 



documented in studies throughout Europe [4-7]. Other 
studies in the US [8-10] or New-Zeeland [11] have shown 
that mortality among COPD patients is directly related 
to a particulate rise in the atmosphere. Here we present 
results concerning elderly males since this age group is 
considered particularly vulnerable to COPD pathology 
[7,12,13]. The link between urban air pollution and COPD 
is well established in time: authors look at daily relations 
between urban particulate concentration and COPD mor- 
tality or morbidity. However there are no analysis devel- 
oping a comparative approach of differences in COPD 
morbidity among places, such as cities. 

Background 

The spatial approach in health research 

There is some methodological debate about how to ap- 
proach the relationship between place and health [14-16]. 
Some authors have pointed out the lack of a firm theoret- 
ical framework for understanding the processes underpin- 
ning this relationship [17,18]. 

In epidemiological studies, space can be regarded in 
two different manners. First, it can be considered as one 
of many potential factors explaining health differences. 
In this case, space refers to the socioeconomic and phys- 
ical characteristics of the area in which the individuals 
studied are located. Different approaches can coexist and 
those that are the most commonly used are based on 
a multilevel modeling. In the second case, the spatial 
health inequalities themselves constitute the main sub- 
ject of the research. The aim is then to understand how 
other geographical patterns are related to spatial health 
differentiations. 

In most multilevel studies, the researchers investigate 
the effects of a place on an individuals health, making a 
distinction between compositional and contextual effects. 
According to Macintyre et al. [17], there are compos- 
itional effects when health differences between places are 
attributable to the under- or over- concentration of indi- 
viduals (classified according to age, social status, type of 
employment, etc.) with a specific health profile. There is 
less consensus on the terminology pertaining to con- 
textual effects. The term context is generally related to 
amenities and opportunity structures in an individuals en- 
vironment. It concerns both the local physical characteris- 
tics in the area of residence (air or water quality, features 
of public spaces such as the presence of parks and road 
networks) and the extent of availability of public or private 
services in these areas, or again the sociocultural and his- 
torical features of the local communities, such as social 
cohesion [2,17]. 

In ecological studies, the focus is on spatial aggregates. 
As Curtis and Jones [16] have pointed out, the individual 
is not always the most relevant unit of analysis, and the 
small-area or regional level are more pertinent for some 



issues. Three types of ecological studies can be distin- 
guished according to the level of observation: 

Classic ecological studies, which investigate spatial 
differences in people's health at a single level of observation 

Often, studies of this type are based on the hypothesis that 
deprivation differences between areas account for a large 
part of spatial health variations. For example, Ben-Shlomo 
et al. [19] demonstrated that male and female mortality 
rates before age 65 in England were strongly associated 
with deprivation (Townsend index) and deprivation in- 
equalities on ward and local authority level. Other ex- 
planatory factors concern the physical characteristics of 
the observed area. Corburn et al. [20], for instance, focus- 
ing on small areas of New-York, examined the effects of 
the built-up environment, housing quality, air pollution 
sources and noxious land use on intra-urban childhood 
asthma hospitalization rates. 

Studies that take into account a geographical level higher 
than that of the object of study 

Indeed, issues relating to compositional and contextual ef- 
fects remain relevant when one considers the meso-level of 
spatial entities: a larger geographical entity, within which 
the observed spatial entity is located, can be regarded as its 
context. For example, Congdon et al. [21] studied health 
and mortality differences between wards in England and 
Wales. They considered the districts as spatial contexts for 
the wards and showed there were significant contextual ef- 
fects for explaining the differences observed between wards 
with similar socioeconomic profiles. In the same manner, 
Phillimore and Morris [22] showed that small areas with 
similar deprivation levels had different mortality rates de- 
pending on their regional location, when the regions 
concerned exhibited different overall levels of deprivation. 
Their hypothesis was that differences in the provision and 
use of health services, in the levels of atmospheric pollut- 
ants and in the types of built-up environment could ex- 
plain these health differences. Although they did not 
explicitly use the term "context", they did refer to differen- 
tiation referring to a higher geographical level, which in 
many studies is considered as a contextual aspect. 

Studies that take into account a geographical level lower 
than that of the object of study 

A number of authors have used an ecological approach to 
consider how far residential segregation can explain differ- 
entiation in health outcomes. Most of such research fo- 
cuses on the relationship between ethnic residential 
segregation and mortality in general, but there are also 
studies on the relationship between socioeconomic resi- 
dential segregation and the incidence of specific illnesses 
such as tuberculosis and heart disease [23-25]. For in- 
stance, in a study involving 47 US cities, Cooper [26] 



Aschan-Leygonie et al. International Journal of Health Geographies 2013, 12:19 
http://www.ij-healthgeographics.eom/content/1 2/1/1 9 



Page 4 of 1 5 



examined cardiovascular mortality among black, Asian 
and white men and women. The author used aggregated 
data at city level Residential segregation in the cities was 
characterized using a Gini coefficient. Cooper [26] showed 
that, independently from mean income variations across 
cities, intra-urban racial residential segregation was asso- 
ciated with higher rates of death from cardiovascular 
disease. Some authors, for example Acevedo- Garcia or 
Subramanian [27,28], also interpret residential segregation 
as a context that affects the health of individuals. 

Aims and hypotheses 

The aim of the present study was to investigate inter- 
urban differences in respiratory health status across the 
55 larger French cities (more than 100,000 inhabitants), 
with particular emphasis on two related issues: 

- What combination (if any) of socioeconomic and 
physical factors is associated with these differences? 

- Is a city's health status related to both its intra- 
urban spatial organization (with possible 
segregation) and its regional context? 

Our first assumption is that a city is a geographical ob- 
ject whose physical and socioeconomic characteristics 
interact with each other to give rise to a specific respira- 
tory health status of the city as an entity. The second 



assumption is that a city's intra-urban organization, i.e. 
spatial organization patterns, likewise contributes to 
forming the respiratory health status of the city as a 
whole. The third assumption is that a city's location in a 
regional context also contributes to its respiratory health 
status. Therefore the physical and socioeconomic charac- 
teristics need to be considered on different geographical 
scales (Figure 1): 

- that of the city with, on the one hand, its 
socioeconomic features (social groups, educational 
levels, health amenities) and its physical features 
(local climate, air pollution, etc.), and on the other 
hand, its intra-urban organizational features, both 
socioeconomic and physical. The hypothesis to be 
examined is for example that two cities with similar 
socioeconomic and physical characteristics but 
different degrees and forms of intra-urban 
residential organization (Figure 2) might differ in 
terms of respiratory health status. 

- that of the region, where physical and 
socioeconomic profiles can be considered as proxies 
for an overall regional context. The hypothesis to be 
explored is that two cities with the same 
socioeconomic and physical characteristics but 
situated in different regional contexts may not have 
the same respiratory health profiles. 



Dunkerque 
Calais,*-^ _j_jlle 

Valenciennes 




Figure 1 Largest French urban units: more than 100,000 inhabitants (2008). 



Aschan-Leygonie et al. International Journal of Health Geographies 2013, 12:19 
http://www.ij-healthgeographics.eom/content/1 2/1/1 9 



Page 5 of 15 




Figure 2 Scales of observation: the regional and city scale. 



The objective of this paper is to analyze the relationship 
between the urban health status, measured for a particular 
pathology and for each city, and its physical and socioeco- 
nomic characteristics, taking into account the regional 
context and the intra-urban organization (Figure 3) b . 

Methods 

In order to work with comparable urban entities, we chose 
to use urban units to define cities. Indeed three deli- 
mitations can be used in France in comparative urban 
studies [29]: the urban central municipality, the urban unit 0 
and the functional urban area d . The urban central munici- 
pality, corresponding to a political grid, is not suited to epi- 
demiological questions. Urban units and functional urban 
areas both offer relevant delimitations for socioeconomic 
comparison. The functional urban area, encompassing ru- 
ral, suburban and urban types of space, covers very differ- 
ent types of physical environment and is therefore not 
suitable. The urban unit, based on morphological continu- 
ity (less than 200 meters separating buildings), covers ex- 
clusively urban built-up areas and is well suited to a 
comparative perspective. As air pollution monitoring is ob- 
ligatory for urban units with more than 100,000 inhabi- 
tants, the 55 French urban units exceeding this threshold 
were included in the study. 

In order to apprehend the intra-urban organization, we 
used the IRIS census subdivisions. These census tracts are 
statistical units defined by INSEE, and known as IRIS 
(Hots Regroupes pour des Indicateurs Statistiques). They 
comprise between 1,800 and 5,000 inhabitants. It is com- 
mon to use this subdivision in France in segregation and 
socio-spatial urban studies [30,31]. 



Data collection 

Literature on the geography of respiratory health has been 
looking at respiratory mortality or hospitalization for 
causes such as lung cancer, COPD or asthma [32-39]. Au- 
thors investigating the link between physical characteris- 
tics, and more specifically pollution and health within 
cities either look at COPD or cardiovascular diseases 
[40,41]. Some authors concentrate on respiratory health 
within cities and most often look at asthma and/or COPD 
[4,6,42-46]. We used the incidence of hospitalization of 
elderly males as a proxy for the respiratory health in a city. 
This data used to create proxies for the urban unit health 
status were derived from the PMSF database which 
gathers information on all hospitalizations, comprising the 
patients place of residence (ZIP code), age, gender, admis- 
sion and discharge dates and all diagnoses. COPD is a 
chronic limitation or obstruction of airflow [12,47] and it 
has been documented as a respiratory disease that is exac- 
erbated by atmospheric pollution [5-7,10,48]. It is often 
considered as an adult (over 25) disability [34] and is more 
frequently diagnosed in males than females. Prevalence 
studies usually consider adults over 20 [12,33,47]. Jeannin 
[13] states that the first COPD hospitalization generally oc- 
curs between ages 65 and 69. Halonen et al. 2013 [7] linked 
COPD hospital emergency room visits by age group (chil- 
dren, adults and elderly > 64 years) with air pollution levels 
and showed a clear effect between COPD hospitalizations 
of elderly and particulate air pollution. Their findings also 
suggest that the mechanisms of respiratory effects of air 
pollution differ by age group. Therefore a study con- 
centrating on only one age group can be interesting. 
According to the results of Tissot-Dupont [49] elderly 
people (> 60 years) are particularly vulnerable to respira- 
tory infectious pathologies. For the present study 12 diag- 
nostic codes describing COPD were chosen according to 
the Furham and Delmas [33] methodology. Data pertaining 
to COPD diagnoses for the year 2008 were extracted to 
identify and assess specific respiratory health patterns. In 
an initial stage we concentrated on all males over 65 
(38,323). Several health indicators were tested: hospita- 
lization rates by age group (65-70; 71-75; 75-80; 
80-85; 86-90) and standardized hospitalization ratio 
(based on standardized mortality ratio methodology) 
for all the age groups pre-cited population. The 65-75 age 
group is considered to be one of the most vulnerable [12] 
for COPD. In the PMSI database, the 71-75 age group 




Figure 3 The intra-urban scale: patterns of spatial organization. 
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is more numerous than the 65-70 (8,340 versus 6,463), 
and the hospitalization ratio is higher (19 versus 15 per 
1,000). Inter-urban patterns for elderly male COPD 
hospitalization rates are very similar to those provided 
by other COPD health indicators (age groups rates and 
standardized hospitalization ratio). Therefore elderly 
(age 71-75) males' COPD hospitalization rates was 
chosen as indicator to characterize the respiratory 
health status of an urban unit. The same indicators 
were calculated for all hospitalizations whatever the 
diagnosis, and this overall hospitalization rate among 
elderly males served as a benchmark that enabled us to 
highlight the specific nature of respiratory health pat- 
terns (Additional file 1: Indicators and measurements). 

A review of the literature on health and socioeconomic 
inequalities [35,42,43,46,50-55] pointed in the direction of 
a large number of indicators relating to various socioeco- 
nomic, physical and amenity aspects. For the socioeco- 
nomic dimension, classic indicators such as unemployment 
rates by age group [17,56-59], the proportion of individuals 
with no qualifications or those having attended at least two 
years of higher education [56,59-62], non-taxable house- 
hold income were established for different scales of obser- 
vation (Additional file 1: Indicators and measurements). 
INSEE f socioeconomic databases dating from 2006 were 
used to create this set of indicators. Family physicians and 
pulmonologists indicators 8 were also established to take 
into account the presence of one important type of health 
care amenity [63,64]. 

In order to characterize the intra-urban organization of 
the urban units in terms of residential segregation, three 
complementary indicators [65,66], were established: the 
coefficient of variation, which was used to represent the 
degree of heterogeneity across census tracts h in the urban 
unit; the Gini concentration ratio, assessing the evenness 
of the spatial distribution of a given characteristic, such as 
educational level, across census tracts; and Morans spatial 
autocorrelation coefficient, which significance level values 
were used to measure the degree of similarity across cen- 
sus tracts. Based on literature review [67-69], educational 
characteristics are more appropriate for the French case 
than ethnic variables. 

A large number of studies have confirmed that respira- 
tory diseases are related to the physical characteristics of 
the living area: COPD has multi-factorial etiology, includ- 
ing exogenous factors like air pollution, cold [70] and top- 
ography [71-73]. Little is however known about the long 
term effects of climate, but there is evidence of statistical 
associations between the prevalence of respiratory diseases 
in different locations and various climatic factors or cli- 
mate zones [11,74-77]. 

Recent scientific papers have discussed the link between 
meteorological and climatic factors and COPD prevalence 
[78]. Across France, because of its geographical situation 



and its rather large surface area (674,800 square kilome- 
ters), there is significant spatial variability in climate. Vari- 
ous climatic zones occur across the country (oceanic, 
semi-continental, mountain, mediterranean). It is there- 
fore reasonable to hypothesize that climate can affect 
inter-urban spatial patterns of COPD prevalence. The 
effect could be direct, for example, an effect of air 
temperature or humidity on the reactivity of the respira- 
tory tract. The effect could also be indirect through a 
differential exposure to air pollutants, infections or aero- 
allergens. If this is the case, spatial patterns of hospital ad- 
missions for COPD might be related to climatic zones and 
variability in climatic parameters like temperature or hu- 
midity. The climatic parameters considered in the litera- 
ture vary: annual mean or seasonal mean temperature 
[11,74-77,79], mean temperature for the coldest and hot- 
test months [75], annual variation of mean temperature 
[76,77,79], relative outdoor humidity [74,76,80,81] fre- 
quency of fog [11], and wind force or direction [11,80]. 
This review of the literature on climatic factors and their 
potential effects on respiratory health (most papers focus 
on asthma) led us to choose a wide array of indicators de- 
scribing the urban units. Firstly, each urban unit was 
assigned to a climate zone which constitutes its regional 
climatic context. On urban unit scale, the following pa- 
rameters were chosen: minimum and maximum tempera- 
tures in January and July, annual mean outdoor relative 
humidity minima, average number of foggy days per year 
(visibility less than one kilometer), average number of days 
with strong winds per year (wind > 57 km/h) and average 
number of hot days per year (> 25°C) were selected. The 
last parameter was chosen because there is clear evidence 
in the literature reporting temporal studies that high 
temperature events (heat waves) increase morbidity in 
general, and in particular for elderly people [82,83]. All 
these variables were derived from Meteo-France meteoro- 
logical data for the period of 1981-2010. 

Two indicators defining the altitude of each urban unit 
were also included: firstly the altitude of the urban unit 
center. There is evidence that COPD prevalence is linked 
to altitude [71-73,84], but there are diverging results. 
Some studies suggest that higher altitude is associated 
with higher COPD prevalence [71,72,84], whereas others 
suggest the reverse relationship [73]. Secondly, in addition 
to the altitude of the urban unit centre, a measure of 
intra-urban variation in altitude was also introduced. 
There is no direct effect of altitude variability on COPD 
prevalence related in the scientific literature. This indica- 
tor was included because it might have an indirect link 
with COPD prevalence via effects on climatic parameters. 

Some studies have investigated associations between 
daily variations in airborne pollen concentrations and re- 
spiratory morbidity or mortality [6,85-87]. These results in- 
dicate a positive relationship between daily rates of COPD 
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mortality and airborne pollen concentrations. Their find- 
ings suggest that other small particles of biological origin 
can potentially have inflammatory effects and exacerbate 
COPD symptoms. In order to take account of the aerobio- 
logical characteristics of the urban units, RNSA 1 pollen risk 
indexes were included. For each pollen type, RNSA defines 
a level of allergic risk and its geographical extension. The 
mapping of the situation in 2008 was used to construct a 
global index, derived from aeroallergen species^ risk levels. 
This index was established for each urban unit. Only pol- 
lens classified as medium, high or very high allergen risk 
by RNSA in 2008 were considered. 

According to the literature there is strong evidence that 
pollutants exacerbate COPD symptoms [4,6,7,39,88-92]. 
Indicators of air pollution were created on the scale of the 
urban unit. Two ambient air pollutants, N0 2 , and PM 10 
were used. The dataset used for assessing the air pollutant 
levels for each pollutant in each urban unit is derived from 
ADEME k and the Geovariances interpolation model. The 
ADEME-Geovariances model was developed from con- 
centrations measured in urban background stations and 
estimated emission data [93]. The model estimates a range 
of pollutant levels: N0 2 (annual mean of daily concentra- 
tions; annual mean of daily 95th percentile concentrations; 
mean of daily winter concentrations) and PM 10 (annual 
mean of daily concentrations; annual mean of the 95th 
percentile of daily concentrations) within cells of a 4x4 
kilometer grid. For each urban unit we calculated a range 
of concentration parameters for the two pollutants, aver- 
aging all cells for each urban unit. As there is a consider- 
able temporal variability in pollutant levels throughout the 
day and the seasons, in addition to the annual mean we 
used parameters that illustrate this variety: annual mean 
of daily 95 th percentile and mean of daily winter (N0 2 ). 
The intra-urban level of variability of air pollutant was 
assessed by calculating a coefficient of variation for each 
urban unit and each pollutant. These variation coefficients 
were calculated from estimated air pollution concentra- 
tions within cells in a grid of 4 km x 4 km. These hetero- 
geneity indicators contribute to the description of citiess 
physical "profile". 

Statistical processing 

Thus a substantial amount of data was mobilized and we 
analyzed it in an exploratory way. To start the relation- 
ships between hospitalization s rates on the one hand, and 
all socioeconomic and physical indicators on the other 
hand, were explored using bivariate statistical analyses 
(correlations and analysis of variance). Secondly, the pur- 
pose was to highlight the complementary roles of different 
indicators liable to account for health status differences 
between urban units, and to assess the impact of each of 
these indicators, ceteris paribus. We proceeded in differ- 
ent steps. First we developed multiple regression analyses 



separately for each family of indicators. We considered 
four families of socioeconomic indicators: unemployment, 
education, medical amenity, the size of the urban unit; 
and three families of physical indicators: climatic and 
topographic features, pollen and finally ambient air pollu- 
tion. This first step was developed in order to select the 
most appropriate indicators in each family and avoid co- 
linearity problems [94] . Indeed, using measurements rela- 
tive to different sub-sets of population (for example 
unemployment for different age groups) on the one hand 
and to different geographical scales on the other (for ex- 
ample unemployment rate both at the urban scale and at 
the regional scale), the indicators were collinear. For each 
family of indicators we selected the ones producing the 
highest explained variance, eliminating those with co- 
linearity. The second step consisted in a forward stepwise 
regression (SAS ©) on all selected indicators. Only classic 
regressions were developed. Indeed, the number of urban 
units by region being very small (21 regions for 55 urban 
units), multi-level analyses were not appropriate. Novelty 
is in the chosen scale of analysis (urban units) and the 
combination of explanatory variables derived from differ- 
ent thematic domains and geographical scales, in an ex- 
ploratory approach. 

Results 

COPD hospitalization rates among elderly men varied by 
a factor of eight, from 6 (Chambery) to 48 per 1,000 indi- 
viduals (Valenciennes). A spatially structured distribution 
was observed, with higher hospitalization rates mostly in 
the urban units in Northern and North-Eastern France, 
and generally lower rates across the rest of France, espe- 
cially in the medium-sized urban units (Figure 4). The 
overall hospitalization rate for all causes ranged from 471 
(Annecy) to 1,162* per 1,000 individuals (Saint-Nazaire), 
with a spatial distribution that was significantly different 
from that of COPD hospitalizations. Most notably, the 
Eastern urban units, some smaller Western ones and 
Marseille indicated high levels, whereas Northern urban 
units showed a combination of medium and low levels. 
The largest urban units showed medium to high levels. 
It appeared clearly that COPD and all hospitalizations 
show different geographies. 

The first step consisted in simple bivariate statistical ana- 
lyses. They showed that most socioeconomic indicators 
were highly correlated with COPD hospitalization rates in 
the urban units (Table 1). As an example, COPD hospitali- 
zations ratio for elderly males were positively and quite 
strongly correlated (r between 0,3 and 0,5) with the per- 
centage of unemployment, the percentage of persons with 
no diplomas, and the percentage of non taxable house- 
holds. On the opposite COPD hospitalizations showed a 
negative correlation with the percentage of individuals with 
high diplomas (r=-0,5). Positive associations (r between 
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COPD hospitalization rates 
among men 71-75 years old 




All causes hospitalization rates 
among men 71-75 years old 



427,659 
227,771 



Administrative regions 



Sources : PMSI 2008, INSEE RGP 2006 
Done with © Philcarto 




] 711 
^ 470.5 



Figure 4 COPD hospitalization rates and overall hospitalization rates among men aged between 71 and 75 in the largest urban units 
in France (2008). 



0,3 and 0,5) were found with the intra-urban organization 
indicators such as the coefficient of variation (CV), or the 
Gini concentration ratio of population with a bachelor or 
bachelor and two years of university college. Negative cor- 
relations (r between -0,3 and -0,5) were found between 
COPD hospitalization rates and intra-urban organization 
indicators such as CV or Gini concentration ratio of popu- 
lation with no diploma. No associations were found with 
the importance of medical amenities or with the size of the 
urban unit (total number of inhabitants). 

Associations between COPD hospitalizations ratio for 
elderly males and physical indicators were much less 
systematic. Only few were significant, and less than the 
correlations between COPD hospitalizations and the socio- 
economic indicators. There was for example a negative as- 
sociation between COPD hospitalizations with January 
maximum temperatures (r=-0,28) and positive with hu- 
midity (r=0,34). There were no associations with other 
physical indicators as altitude, pollen or wind. Regarding 
air pollution, there were a weak positive correlation (r 
around 0,3) with PM 10 annual mean of daily concentra- 
tions and annual mean of the 95th percentile of daily con- 
centrations. A nonlinear relationships was observed 
between N0 2 levels and COPD hospitalization rates. In 
most urban units, respiratory health status seemed to de- 
teriorate as the N0 2 level increased. However, the highest 
N0 2 concentrations correspond to urban units with mod- 
erate COPD hospitalization rates, and not with the highest, 
as could have been expected (Figure 5). 

Conversely, there was almost no relationship between 
the inter-urban distribution of the benchmark indicator 
"overall hospitalizations" and the socioeconomic and 



physical indicators. The only correlations that were statis- 
tically significant concerned the CV and Gini coefficient 
defining the intra-urban distribution of populations 
according to their levels of education. These correlations 
were positives, and thus similar to those observed for 
COPD hospitalizations. Except for these indicators, nei- 
ther socioeconomic nor physical indicators were corre- 
lated with "overall hospitalizations". 

In a second step we developed multiple regressions sep- 
arately for each family of indicators in order to select the 
most appropriate for the final analysis. In some cases only 
one indicator was selected. This was the case for un- 
employment. The variable producing the highest ex- 
plained variance was the rate for 15-64 unemployment 
age group at the regional scale. Once this indicator taken 
into account, no other indicator belonging to the un- 
employment family of indicators was significant. In other 
cases a few indicators were selected. As an example, we 
found an additive effect of the two pollutants that were in- 
vestigated, whereby together they accounted for 33% of 
the inter-urban variations in COPD rates. 

In a third step a forward stepwise regression was 
performed on the selected indicators in preceding step in 
order to explore the relationship between respiratory 
health status and the different indicators studied ceteris 
paribus. The analyses revealed a combination of factors 
that explained the statistics pertaining to health status in 
the urban units. The model obtained for elderly male com- 
bined socioeconomic and physical indicators on different 
geographical scales: an unemployment indicator on re- 
gional scale, a climate indicator on urban scale and two in- 
dicators reflecting intra-urban residential organization. 
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Table 1 Relationships between hospitalizations and some of the explanatory indicators for COPD and overall 
hospitalizations 

Type of indicator COPD hospitalizations Overall hospitalizations 



Signi- Type of relationship Signi- Type of relationship 

ficance ficance 
level level 



Socio- 


Urban unit scale : 


n.s 


no significant relation 


n.s 


no significant relation 


economic 
indicators 


Size (number of inhabitants) 












uruun uf lit bLUic . 


n nm 

U.UU I 


LI Ileal. pilVlleLjeU LllleS Nave lUVVel 


n.s 


IIU blLjl III ILdl IL leldllUII 




Socioeconomic (% unemployment/ft of 




hospitalization rates 








individuals with no diploma, bachelor 




Positive with unemployment, non 








diploma,% of non taxable 




taxable households (0.3<r<0.5) 








households, etc.) 




Negative with bachelor diploma 
(r=-0.5) 








Intra-urban scale 


0.001 


Linear positive: cities with 


0.001 


Linear positive: cities with 




Residential intra-urban organization 
indicators (CV, Gini, Moran of 
unemployment/^ of individuals with no 




concentration of privileged 
populations have higher 
hospitalization rates (0.3<r<0.44) 




concentration of privileged 
populations have higher rates 
(0.28<r<0.40) 




diploma, bachelor diploma,% of non 




Linear negative: cities with 




Linear negative: cities with 




taxable households, etc.) 




concentration of deprived populations 
have lower rates (-0.3<r<-0.48) 




concentration of deprived 
populations have lower rates 
(r=-0.43) 




Regional scale 


0.001 


Linear positive: cities in deprived 


n.s 


no significant relation 




Socioeconomic (% unemployment, etc.) 




regions have higher hospitalization 
rates (0.33<r<0.47) 








Urban unit scale 


n.s 


no significant relation 


n.s 


no significant relation 




Access to health amenities 










Physical 


Urban unit scale : 


0.05 


Linear negative with the daily mean 


0.05 


no significant relation 


indicators 


Climate parameters (temperature, 
humidity) 




temperature in January (r=-0.28) 

Linear positive with the mean annual 
minimal relative humidity (r=0.34) No 
significant relation with other climate 
indicators 








Regional scale 


n.s 


no significant relation 


0.007 


Higher hospitalization rates in 




Climate zones 








semi-continental climate ; 
lower rates in mountain 
climate 




Urban unit scale : 


0.05 


Globally lower hospitalization rates for 


n.s 


no significant relation 




Air pollution (ambient levels of PM 10 , 




less polluted locations 








N0 2 ) 




Linear positive for PM 10 (0.28<r<0.30 ) 
Non linear positive for N0 2 








Intra-urban scale : 


n.s 


no significant relation 


n.s 


no significant relation 




Air pollution (intra-urban variations of 












PM 10 , N0 2 ) 












Urban unit scale 


n.s 


no significant relation 


n.s 


no significant relation 




Pollen index 











This model explained 53% of inter-urban differences in 
COPD hospitalization rates (Table 2). This R 2 was high 
compared to the 27% obtained with the model derived for 
"overall hospitalizations", where only two indicators (both 
relating to the intra-urban socioeconomic organization) 
were significant. A comparison of the results from these 
models suggests that socioeconomic and physical factors 
are both specifically discriminating for respiratory health 
status. The model residuals showed no spatial auto- 



correlation. We verified the robustness of these results on 
the standardized indexes. 

Discussion 

Lower levels for COPD hospitalization ratios for elderly 
males tended to be found in cities located in Western and 
South- Western regions. High levels were mostly found in 
Northen and North-Eastern regions. This spatial pattern 
is consistent with that obtained by Delmas and Fuhrman 
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Figure 5 Nonlinear relationship between average COPD hospitalization rates among elderly men and nitrogen dioxide concentration 
levels (2008). 



in 2010 [32] on COPD and asthma at the regional scale 
and by Rican et al [36] in their 2003 study on respiratory 
diseases in French urban units with more than 20,000 in- 
habitants. These authors suggest that there might be an 
effect of mining and steel manufacturing industry. This is 
true for several urban units of this study (St-Etienne, 
Thionville and Bethune). However other urban units like 
Creil, Caen or Calais that have no such specialization also 



show high rates. This issue needs to be further 
investigated. 

The geography of both COPD and overall hospitaliza- 
tions is not specific to any age group, since the stan- 
dardized indexes show the same patterns. The proxies 
used to describe urban unit health status could be im- 
proved. Indeed, in this study, they were based on hospi- 
talizations alone, which, according to Vigneron [95], 



Table 2 Results for the multiple regression models explaining respectively for "COPD hospitalization" rates and 
"overall hospitalization" rates 



Explained variable 


Explanatory variables 


Coefficient 


Significance level 


COPD hospitalization rates for men aged 
between 71 and 75, calculated as proportion 
of male population in same age group 


Intercept 

Unemployment rates for 15-64 age group 
(Regional scale) 


11 
171 


0.24 
0.002 




Gini coefficient (Intra-urban scale): 

individuals with bachelor degree and two years 

of college university study as proportion of overall 

population 


54 


0.023 




Coefficient of variation (Intra-urban scale) 
of individuals with no diploma 


-35 


0.0006 




Daily temperature maximum in January 
(Urban unit scale) 


-0,93 


0.037 




% of explained variance (R 2 ) 


53% 


0.0001 


Explained variable 


Explanatory variables 


Coefficient 


Significance level 


Overall hospitalization rates for men aged 
between 71 and 75, calculated as proportion 
of male population in same age group 


Intercept 

Gini coefficient (Intra-urban scale): 

individuals with bachelor degree and two years 

of college university study as proportion of overall 

population 


1002 
911 


0.0001 
0.04 




Coefficient of variation (census tracts scale) 
of non-taxable households 


-866 


0.001 




% of explained variance (R 2 ) 


27% 


0.002 
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only represents half of total health care consumption in 
France. Tonnellier [96] showed that hospital consulta- 
tions are more frequent than visits to the family phys- 
ician in Northern France, as well as in medium-sized 
and smaller urban units. We intend to conduct further 
investigations to gain better insight into the inter-urban 
differences in respiratory health status, by considering 
both hospitalizations and visits to family physicians and 
specialists. 

Multiple regression analyses were used to explain inter- 
urban differences for both COPD and overall hospi- 
talization rates. Many statistically significant relationships 
were observed between COPD rates and the candidate ex- 
planatory socioeconomic and physical dimensions. The 
COPD models combined indicators associated with both 
physical and socioeconomic characteristics and involved 
different geographical scales. Conversely, very few strong 
relationships were observed for overall hospitalization 
rates. This reinforces our hypothesis that respiratory 
health status exhibits specific inter-urban patterns. 

Our first hypothesis that the socioeconomic and physical 
dimensions both need to be considered in order to provide 
a statistical explanation for COPD differences across urban 
units is not clearly borne out. Our results showed that 
inter-urban differences in respiratory health status are 
more systematically related to socioeconomic factors than 
physical ones. The relationships between health and socio- 
economic indicators were similar to those established using 
a deprivation index in intra-urban studies [36,53,97,98]. 
These results are also consistent with studies on socioeco- 
nomic inequalities and health carried out at different geo- 
graphical scales, for example by Leclerc et al. [53] at a 
regional scale, Rican, et al. [35,36] at the city and regional 
scale, Blomgren et al. [99] at the regional scale, and 
Macintyre et al. [98] at the small-area (intra-urban) scale. 

The positive correlation between average annual humid- 
ity and COPD hospitalization rates confirmed the evidence 
reported in the literature [11,76,80,81]. One of the climate 
parameters integrated into the model contributed clearly 
to explaining inter-urban differences: the annual mean of 
maximum temperature in January was negatively associ- 
ated with hospital admission rates for COPD. This is 
consistent with results that have shown the existence 
of a relationship between temperatures and asthma 
[11,75-77,79]. However, for asthma prevalence the direc- 
tion of the relationship is positive with annual temperature 
means [11,77,79] or lowest monthly temperature means in 
the ISAAC study [76]. Our results are different, as they 
showed that the higher the temperature in January, the 
lower the COPD rates. A direct comparison of results is 
not possible, because in the literature asthma was the ob- 
ject of study and not COPD prevalence, and the age groups 
observed were different: children [76], young adults or 
adults [11,75,77,79]. However, the negative relationship 



between temperature and hospital admissions for COPD 
clearly appears in temporal studies and it is in January that 
the maxima of monthly hospital admissions are observed 
[78]. There is no confirmation in the literature that the as- 
sociation between temperature and respiratory health out- 
comes reflects a causal relationship, or that correlations are 
a result of indirect relationships and linked to other factors, 
like air pollution levels. 

The adverse effect of air pollution on health has been 
well established in literature, and our results are consistent 
with only some of this evidence. COPD hospitalization 
rates were significantly, but weakly correlated, with all 
parameters describing PM 10 concentrations in the ur- 
ban units. This is consistent with the literature, where 
COPD appears clearly linked with PM 10 ambient levels 
[7,39,88-92]. The non linearity of the relation between 
N0 2 levels and COPD hospitalization rates of elderly 
males could be attributed to a compensatory effect 
brought about by the higher socioeconomic status of these 
larger urban units, which are also characterized by the 
highest air pollution levels. This is consistent with the 
conclusions of Cakmak et al. [100] in their study on the 
way community income and education modify the effect 
of gaseous air pollution on respiratory hospitalizations. In- 
deed, based on analyses of daily time series for respiratory 
hospitalizations and daily concentrations of N0 2 concen- 
trations in 10 large Canadian cities, they found that the ef- 
fect of both N0 2 and a combination of pollutants was 
stronger with decreasing levels of household income. 

However, in our regression model the air pollutant indi- 
cators did not appear. In our study, socioeconomic indica- 
tors seem to obliterate the effect of nitrogen dioxide 
variations on COPD hospitalizations. This partly invali- 
dates our earlier hypothesis regarding the contribution of 
atmospheric pollution to the inter-urban differences in 
COPD rates: once socioeconomic and climate indicators 
are taken into account, the concentration of pollutants 
such as nitrogen dioxide and particulate matter no longer 
provides a statistical explanation for variations in COPD 
rates across urban units. 

Our second hypothesis, that explanatory factors should 
be considered on different spatial scales to account for dif- 
ferences in respiratory health status among urban units, 
was confirmed. The most innovative results were those 
that involved the use of regional and intra-urban indica- 
tors. First, when combined with other variables, un- 
employment showed a more significant contribution on 
regional than on urban scale. Thus, an urban units re- 
spiratory health status could be partially explained by its 
belonging to a regional economic context. This result is 
consistent with the findings of Rican et al. [36,53,97,98], 
who reported that the regional component of respiratory 
mortality variations in France is highly significant. It is 
also consistent with Jusots result [101] that the individual 
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mortality risk among male adults in France includes a re- 
gional component. Secondly, an intra-urban organization 
effect was clearly observed. In the literature, it has been 
shown that there are positive associations between marked 
racial residential segregation and mortality rates or to- 
bacco risk on city scale [23-25] . Several explanations have 
been put forward: segregation may limit the social and 
economic opportunities for deprived groups, worsen their 
deprivation, and in turn lower the city s overall health sta- 
tus; similarly, residential segregation could restrict mobil- 
ity and access to health services, and consequently affect a 
city s health status. We expected to observe the same over- 
all effect of marked residential differentiation on health 
status. However, the statistical results proved far from 
simple: residential organization of the undereducated 
population within the urban unit was found to correlate 
negatively with COPD rates on urban unit scale; con- 
versely, the residential organization of highly educated 
populations was found to correlate positively with COPD 
rates. This seems to imply that the spatial concentration 
of less educated populations is associated with better re- 
spiratory health status, whereas spatial concentration of 
more highly educated populations is associated with lower 
respiratory health status at the urban unit scale. Even 
though the number of IRIS census tracts varied between 
urban units, no statistical effect arising from differences in 
urban unit sizes seems to have affected these results. The 
relationships with residential organization are then par- 
tially coherent with the conclusions of various studies 
cited above, in the sense that the results are far from sim- 
ple and will require further investigation. 

Conclusion 

It is important to gain a better understanding of the 
spatial patterns of urban respiratory health problems, 
since respiratory health is a growing cause of morbidity 
worldwide. It is reported as being dependent on both so- 
cioeconomic and physical environmental factors, more 
specifically pollutants. Our aim is to explore these hypoth- 
eses at the inter-urban scale. We do not aim to provide 
any predictive model about a disease outcome and envir- 
onmental exposure, neither to establish any causality, at 
the level of the city or at the individual level. 

First the study shows surprising results for the contribu- 
tion of pollution factors. Indeed, differences in respiratory 
health among cities related largely to socioeconomic fac- 
tors and only marginally to pollution. Second, the study 
underlines that broader contextual effects come into play. 
Indeed it appeared that differences in health respiratory 
status between cities arose from the economic regional 
context (in terms of unemployment ): thus characteristics 
observed at the regional scale contribute to explain inter- 
urban differences in respiratory health status. Third, 
inter-urban inequalities in health status are related to 



intra-urban organizational factors, such as differentiated 
residential locations of specific population groups (highly 
qualified groups or non-qualified groups). In order for pub- 
lic policy makers to plan "healthy cities" and thus improve 
overall public health, it is essential to gain a better under- 
standing of the complex and probably very interdependent 
factors that contribute to the development of spatial health 
variations. Indeed, the hypotheses examined - 1) that two 
cities with similar socioeconomic and physical environ- 
ment characteristics but different degrees of intra-urban 
residential organization could differ in terms of respiratory 
health status and 2) that two cities with the same socioeco- 
nomic and physical characteristics but situated in different 
regional contexts may not have the same respiratory health 
profiles - are confirmed. 

Endnotes 

a Our aim is by no means to state that an inhabitant in 
city X has a greater or lesser risk of contracting a respira- 
tory disease than an inhabitant in city Y. 

b Each shade of grey represents the density of presence 
of a certain type of population (highly educated persons 
for example): dark-grey stands for a high density, pale- 
grey for a low density. The examples depict common 
intra-urban patterns of residential organization. These dif- 
ferences form segregation patterns within the city: from 
center towards periphery in the first graph, marked con- 
trast in the second (in these two graphs intra-urban differ- 
entiation is high), lesser contrast in the third and fourth 
graph (intra-urban differentiation is lower). 

c An Urban Unit is a municipality or group of municipal- 
ities with at least 2,000 inhabitants. An urban unit is a con- 
tinuous built up area in the sense that there is less than 
200 meters between two buildings (Insee web site, 2013). 

d The Functional Urban Area can be defined as travel- 
to-work area. Principally it is an agglomeration of work 
places attracting the work force from the surrounding area. 

e PMSI: Programme Medical Systeme dlnformation 
(Source: Agence Technique de Tlnformation Hospitaliere 
[http://www.atih.sante.fr]). 

f INSEE ("Institut National de la Statistique et des 
Etudes Economiques") is the French National Institute 
of Statistics and Economic Studies (http://www.insee.fr). 

g Data is derived from the "Base Permanente des 
Equipements" at INSEE and originally from ADELI (Au- 
tomatisation DES LIstes), which keeps records of all regis- 
tered health practitioners in France. 

h These census tracts are statistical units defined by 
INSEE, known as IRIS (Hots Regroupes pour des Indi- 
cateurs Statistiques). They roughly comprise between 
1,800 and 5,000 individuals. 

1 "Reseau National de Surveillance Aerobiologique". 

] The ten pollen species used to build the aerobiological 
indicator are : 
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- Trees: alder (Alnus), cypress (Cupressus), oak 
{Quercus), ash (Fraxinus), birch (Betula), plane 
{Platanus), olive (Olea) 

- Grasses {Family Poaceae): no distinction was made 
between species of grass 

- Weeds: ragweed {Ambrosia), nettle/parietaria 
{Urticaceae) 

k ADEME, "Agence de TEnvironnement et de la 
Maitrise de TEnergie", is the French Environment and En- 
ergy Management Agency [www2.ademe.fr]. 

1 There can be more than 1,000 hospitalization for 1,000 
individuals as one individual can have paid more than one 
visit to the hospital in a year. 
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